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Avian coronavirus infectious bronchitis virus (IBV) is variable, which causes many serotypes. Here we reported the complete 
genome sequences of two virulent IBV variants from China, GX-YL5 and GX-YL9, belonging to different serotypes. Differences 
between GX-YL5 and GX-YL9 were found mainly in stem-loop structure I in the predicted RNA secondary structure of open 
reading frame (ORF) 1b and the S protein gene fusion region, which will help us understand the molecular evolutionary mecha- 
nism of IBV and the disconcordance between the genotypes and serotypes of coronavirus. 


‘aie bronchitis virus (IBV), a member of the genus Gam- 
macoronavirus of the family Coronaviridae in the order Nidovi- 
rales, is an enveloped, nonsegmented, single-stranded, positive- 
sense RNA virus. The coronavirus employs a discontinuous RNA 
synthesis mechanism for its transcription, and multiple sub- 
genomic mRNAs contain a short 5’ leader sequence. The leader 
and mRNA body sequences are joined within a short conserved 
sequence motif named the transcription-regulating sequence, 
which precedes each transcription unit (4, 8, 10, 12). In this study, 
two IBV field strains from China, named GX-YL5 and GX-YL9, 
showed strong pathogenicity and belonged to the same genotype 
but different serotypes (5). The complete genomes of GX-YL5 and 
GX-YL9 were sequenced and analyzed to further understand the 
molecular mechanism of IBV. 

The 5’ and 3’ terminals of the GX-YL5 and GX-YL9 genomes 
were confirmed by using the 5’-Full random amplification of 
cDNA ends (RACE) kit and 3’-Full RACE core set, version 2.0 
(TaKaRa, Japan), and the other parts were generated by 25 and 24 
overlapping cDNA fragments, respectively. All sequencing was 
done using an ABI Prism 3730 sequencer (Applied Biosystems) 
and assembled using the SeqMan software program (DNASTAR 
Inc.). Sequence alignment was conducted and a phylogenetic tree 
was constructed using the software program MEGA5 (3). RNA 
secondary structure prediction was conducted using the program 
Mfold (15). Recombination analysis was performed using the 
RDP 4.14 (7) and SimPlot 3.5.1 (6) software programs. 

The complete genomes of GX-YL5 and GX-YL9 are 27,706 and 
27,582 nucleotides (nt), respectively, in length, after excluding the 
polyadenylated tract, which shared 95.9% nucleotide identity. 
Compared with the most popularly used vaccine strain, H120, 
they had lower nucleotide identities, 86.0% and 85.9%, respec- 
tively. Their genome organizations are classical IBV genomes with 
the characteristic gene order 5’-Pol-S-3a-3b-E-M-5a-5b-N-3’. 
The differences between GX-YL5 and GX-YL9 were found in the 
S1 subunit within nt 20429 to 20830, a membrane gene within nt 
24559 to 25579, and the 3’ untranslated region (UTR) within nt 
27193 to 27405. The recombination sites of genes $2, membrane, 
nucleocapsid, and 3’ UTR in GX-YL5 and the recombination sites 
of genes 3a, membrane, 5a, nucleocapsid, and 3’UTR in GX-YL9 
had been predicted. 

The S1 subunit of the IBV genome is the major determiner of 
serotype (1, 2, 9, 13). In our study, the predicted RNA secondary 
structure of the ORF 1b and the S protein gene fusion region (nt 
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20375 to 20430), based on GX-YL5, showed that the nucleotide 
mutations were found mostly in stem-loop structure I, and mul- 
tiple nucleotide differences were also found in the fusion region 
between GX-YL5 and GX-YL9. Interestingly, GX-YL5 and GX- 
YL9 belonged to the same genotype but different serotypes (5). 
During synthesis of the cDNA strand of coronavirus, a process 
known as copy choice can lead to genetic recombination (14). The 
high frequency of recombination in IBV likely plays a major role 
in the generation of new serotypes of the virus (11). 

The results will help us understand the molecular evolutionary 
mechanisms of IBV and disconcordance between the genotypes 
and serotypes of coronaviruses. 

Nucleotide sequences accession numbers. The full genomic 
sequences of GX-YL5 and GX-YL9 are available in GenBank un- 
der accession numbers HQ848267 and HQ850618, respectively. 


ACKNOWLEDGMENTS 


This work was supported by grants from the National Natural Science 
Foundation of China (30700599 and 31160516), the Guangxi Natural 
Science Foundation (0991044), and the Guangxi Provincial Programs 
for Science and Technology Development (0993009-2). 


REFERENCES 


1. Cavanagh D, et al. 1992. Location of the amino acid differences in the S1 
spike glycoprotein subunit of closely related serotypes of infectious bron- 
chitis virus. Avian Pathol. 21:33—43. 

2. Cavanagh D, Davis PJ, Mockett AP. 1988. Amino acids within hyper- 
variable region 1 of avian coronavirus IBV (Massachusetts serotype) spike 
glycoprotein are associated with neutralization epitopes. Virus Res. 11: 
141-150. 

3. Cobo EDC, Silveira TP, Micheletti AM, Crema E, Adad SJ. 2012. 
Research on Trypanosoma cruzi and analysis of inflammatory infiltrate in 
esophagus and colon from chronic chagasic patients with and without 
mega. J. Trop. Med. 2012:232646. doi:10.1155/2012/232646. 

4. Lai MM, Cavanagh D. 1997. The molecular biology of coronaviruses. 
Adv. Virus Res. 48:1-100. 

5. LiM, et al. 2012. Serotype and genotype diversity of infectious bronchitis 


Received 20 July 2012 Accepted 20 July 2012 


Address correspondence to Meilan Mo, momeilan@163.com, or Ping Wei, 
pingwei8@126.com. 


M.M. and B.H. contributed equally to this article. 
Copyright © 2012, American Society for Microbiology. All Rights Reserved. 
doi:10.1128/JVI.01895-12 


jviasm.org 10903 


NSGN Aq GLoz ‘92 Ae UO /Bio'wse'lAl//:dyy wos papeojumoq 


Genome Announcement 


viruses isolated during 1985-2008 in Guangxi, China. Arch. Virol. 157: 
467-474. 

. Lole KS, et al. 1999. Full-length human immunodeficiency virus type 1 
genomes from subtype C-infected seroconverters in India, with evidence 
of intersubtype recombination. J. Virol. 73:152-160. 

. Martin D, Rybicki E. 2000. RDP: detection of recombination amongst 

aligned sequences. Bioinformatics 16:562—563. 

. Nakaya N, et al. 2011. Low-dose pravastatin and age-related differences 

in risk factors for cardiovascular disease in hypercholesterolaemic Japa- 

nese: analysis of the management of elevated cholesterol in the primary 

prevention group of adult Japanese (MEGA study). Drugs Aging 28:681- 

692. 

. Niesters HG, Bleumink-Pluym NM, Osterhaus AD, Horzinek MC, van 

der Zeijst BA. 1987. Epitopes on the peplomer protein of infectious 

bronchitis virus strain M41 as defined by monoclonal antibodies. Virology 

161:511-519. 


10904 = jvi.asm.org 


0. 


Snijder EJ, Meulenberg JJ. 1998. The molecular biology of arteriviruses. 
J. Gen. Virol. 79:961-979. 


. Thor SW, Hilt DA, Kissinger JC, Paterson AH, Jackwood MW. 2011. 


Recombination in avian gamma-coronavirus infectious bronchitis virus. 
Viruses 3:1777-1799. 


. van der Most RG, de Groot RJ, Spaan WJ. 1994. Subgenomic RNA synthesis 


directed by a synthetic defective interfering RNA of mouse hepatitis virus: a 
study of coronavirus transcription initiation. J. Virol. 68:3656—-3666. 


. Wang L, Junker D, Hock L, Ebiary E, Collisson EW. 1994. Evolutionary 


implications of genetic variations in the $1 gene of infectious bronchitis 
virus. Virus Res. 34:327-338. 


. Weiss RA, Mason WS, Vogt PK. 1973. Genetic recombinants and 


heterozygotes derived from endogenous and exogenous avian RNA tumor 
viruses. Virology 52:535-552. 


. Zuker M. 2003. Mfold web server for nucleic acid folding and hybridiza- 


tion prediction. Nucleic Acids Res. 31:3406-3415. 


Journal of Virology 


NSGN Aq GLog ‘92 Ae Uo /Bio'wselAl//:dyy wos papeojumoq 


